Agile parallel bioinformatics workflow management using Pwrake

<p>Abstract</p> <p>Background</p> <p>In bioinformatics projects, scientific workflow systems are widely used to manage computational procedures. Full-featured workflow systems have been proposed to fulfil the demand for workflow management. However, such systems tend to be over-weighted for actual bioinformatics practices. We realize that quick deployment of cutting-edge software implementing advanced algorithms and data formats, and continuous adaptation to changes in computational resources and the environment are often prioritized in scientific workflow management. These features have a greater affinity with the agile software development method through iterative development phases after trial and error.</p> <p>Here, we show the application of a scientific workflow system Pwrake to bioinformatics workflows. Pwrake is a parallel workflow extension of Ruby's standard build tool Rake, the flexibility of which has been demonstrated in the astronomy domain. Therefore, we hypothesize that Pwrake also has advantages in actual bioinformatics workflows.</p> <p>Findings</p> <p>We implemented the Pwrake workflows to process next generation sequencing data using the Genomic Analysis Toolkit (GATK) and Dindel. GATK and Dindel workflows are typical examples of sequential and parallel workflows, respectively. We found that in practice, actual scientific workflow development iterates over two phases, the workflow definition phase and the parameter adjustment phase. We introduced separate workflow definitions to help focus on each of the two developmental phases, as well as helper methods to simplify the descriptions. This approach increased iterative development efficiency. Moreover, we implemented combined workflows to demonstrate modularity of the GATK and Dindel workflows.</p> <p>Conclusions</p> <p>Pwrake enables agile management of scientific workflows in the bioinformatics domain. The internal domain specific language design built on Ruby gives the flexibility of rakefiles for writing scientific workflows. Furthermore, readability and maintainability of rakefiles may facilitate sharing workflows among the scientific community. Workflows for GATK and Dindel are available at <url>http://github.com/misshie/Workflows</url>.</p

Copy Number Alteration and Uniparental Disomy Analysis Categorizes Japanese Papillary Thyroid Carcinomas into Distinct Groups

The aim of the present study was to investigate chromosomal aberrations in sporadic Japanese papillary thyroid carcinomas (PTCs), concomitant with the analysis of oncogene mutational status. Twenty-five PTCs (11 with BRAFV600E, 4 with RET/PTC1, and 10 without mutation in HRAS, KRAS, NRAS, BRAF, RET/PTC1, or RET/PTC3) were analyzed using Genome-Wide Human SNP Array 6.0 which allows us to detect copy number alteration (CNA) and uniparental disomy (UPD), also referred to as copy neutral loss of heterozygosity, in a single experiment. The Japanese PTCs showed relatively stable karyotypes. Seven cases (28%) showed CNA(s), and 6 (24%) showed UPD(s). Interestingly, CNA and UPD were rarely overlapped in the same tumor; the only one advanced case showed both CNA and UPD with a highly complex karyotype. Thirteen (52%) showed neither CNA nor UPD. Regarding CNA, deletions tended to be more frequent than amplifications. The most frequent and recurrent region was the deletion in chromosome 22; however, it was found in only 4 cases (16%). The degree of genomic instability did not depend on the oncogene status. However, in oncogene-positive cases (BRAFV600E and RET/PTC1), tumors with CNA/UPD were less frequent (5/15, 33%), whereas tumors with CNA/UPD were more frequent in oncogene-negative cases (7/10, 70%), suggesting that chromosomal aberrations may play a role in the development of PTC, especially in oncogene-negative tumors. These data suggest that Japanese PTCs may be classified into three distinct groups: CNA+, UPD+, and no chromosomal aberrations. BRAFV600E mutational status did not correlate with any parameters of chromosomal defects

Mutations in UVSSA cause UV-sensitive syndrome and impair RNA polymerase IIo processing in transcription-coupled nucleotide-excision repair

UV-sensitive syndrome (UVSS) is a genodermatosis characterized by cutaneous photosensitivity without skin carcinoma1, 2, 3, 4. Despite mild clinical features, cells from individuals with UVSS, like Cockayne syndrome cells, are very UV sensitive and are deficient in transcription-coupled nucleotide-excision repair (TC-NER)2, 4, 5, which removes DNA damage in actively transcribed genes6. Three of the seven known UVSS cases carry mutations in the Cockayne syndrome genes ERCC8 or ERCC6 (also known as CSA and CSB, respectively)7, 8. The remaining four individuals with UVSS, one of whom is described for the first time here, formed a separate UVSS-A complementation group1, 9, 10; however, the responsible gene was unknown. Using exome sequencing11, we determine that mutations in the UVSSA gene (formerly known as KIAA1530) cause UVSS-A. The UVSSA protein interacts with TC-NER machinery and stabilizes the ERCC6 complex; it also facilitates ubiquitination of RNA polymerase IIo stalled at DNA damage sites. Our findings provide mechanistic insights into the processing of stalled RNA polymerase and explain the different clinical features across these TC-NER–deficient disorders

Behavior of geodesic acoustic mode and limit-cycle oscillation approaching L-H transition in JFT-2M tokamak

In this paper, a phenomenology of competing behavior between the geodesic acoustic mode (GAM) and the limit-cycle oscillation (LCO) is presented. Before the LCO occurs, the GAM can grow to the observable amplitude via the turbulent Reynolds stress force. Approaching the L-H transition, the LCO is excited and the GAM decays. In the LCO phase, the GAM driving force is possibly suppressed by the nonlocal turbulence amplitude modulation by the LCO